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ABSTRACT 

X 

We report the detection of numerous far-infrared emission lines of water vapor 
toward the supergiant star VY Canis Majoris. A 29.5 — 45 (j,m. grating scan of VY 
CMa, obtained using the Short Wavelength Spectrometer (SWS) of the Infrared Space 
Observatory (ISO) at a spectral resolving power A/AA of approximately 2000, reveals 
at least 41 spectral features due to water vapor that together radiate a total luminosity 
~ 25 Lq. In addition to pure rotational transitions within the ground vibrational 
state, these features include rotational transitions within the (010) excited vibrational 
state. The spectrum also shows the 2 n 1/ / 2 (J = 5/2) Yl 3 / 2 (J = 3/2) OH feature 
near 34.6 /im in absorption. Additional SWS observations of VY CMa were carried out 
in the instrument's Fabry-Perot mode for three water transitions: the 7 2 5 — 6i6 line at 
29.8367 /im, the 441 — 3i2 line at 31.7721 /iin, and the 432 — 3o3 line at 40.6909 ^m. 
The higher spectral resolving power A/AA of approximately 30,000 thereby obtained 
permits the line profiles to be resolved spectrally for the first time and reveals the "P 
Cygni" profiles that are characteristic of emission from an outflowing envelope. 
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1. Introduction 

A wealth of recent data from the Infrared Space Observatory (ISO; Kessler et al. 1996) has 
demonstrated that water vapor is a ubiquitous constituent of the astrophysical universe. With 
its powerful spectroscopic capability unhindered by atmopheric absorption and covering the 
2.38 — 197 /im spectral region without interruption, ISO has detected water vapor in astrophysical 
environments as diverse as star-forming interstellar clouds (e.g. van Dishoeck & Helmich 1996; 
Helmich et al. 1996; Ceccarelli et al. 1998; van Dishoeck et al. 1998); Herbig Haro ojects (Liseau 
et al. 1996); shocked interstellar regions (e.g. Harwit et al. 1998); the atmospheres of Jupiter, 
Saturn, Uranus, Neptune and Titan (Feuchtgruber et al. 1997; Lellouch et al. 1997; Coustenis et 
al. 1998); the photospheres of cool stars (e.g. Tsuji et al. 1997); and the circumstellar outflows 
from oxygen-rich evolved stars (e.g. Neufeld et al. 1996; Barlow et al. 1996, Justtanont et al. 
1996). 

The data now available from ISO observations support earlier theoretical predictions of a 
substantial abundance of water vapor within warm molecular astrophysical environments where 
the elemental abundance of oxygen exceeds that of carbon. Since the water molecule possesses a 
large dipole moment, water rotational transitions play an important role in the thermal balance of 
the gas (e.g. Neufeld, Lepp & Melnick 1995). 

As part of a program to study water vapor in shocked interstellar regions and circumstellar 
outflows, we have carried out a search for far- infrared rotational emissions from water toward VY 
CMa, a high-mass red supergiant of spectral type M5Iae that is believed to be in an advanced 
stage of evolution with a bolometric luminosity ~ 5 x 10 5 L and a mass loss rate that exceeds 
lO~ 4 M yr _1 (Danchi et al. 1994). A complete spectrum of VYCMa covering the 2.38 to 45 /jm 
spectral region was obtained using the Short Wavelength Spectrometer (SWS; de Grauuw et al. 
1996) on board ISO and will be presented in a future publication; in this Letter we consider only 
the 29.5 - 45 /im region (SWS Band 4) that contains the strongest rotational emissions from water. 

2. Observations 

All the observations of VY CMa reported here were carried out using the SWS instrument. A 
complete spectrum covering the 2.38 to 45 /im spectral region was obtained using the instrument's 
grating mode at a spectral resolving power A/AA ~ 1000 — 2500 (the 29.5 - 45 /jm portion of 
which is presented here.) In addition, the spectra of three rotational lines of ortho-water were 
obtained in Fabry-Perot (FP) mode at a spectral resolving power A/AA ~ 3 x 10 4 : the 7 25 - 6i 6 
line at 29.8367 /mi, the 4 4i - 3 i2 line at 31.7721 /im, and the 4 32 - 3 03 line at 40.6909 /im. 

The observations were all carried out on 1997 November 19 (grating scan) and 1997 November 
20 (FP scans) with the ISO beam centered at position a = 7 h 22 m 58.3 s , 5 = -25°46'03". The beam 
size was 20 x 33" for the grating scan and 17 x 40" for the FP scan. Including overheads for dark 
current measurements and calibration, the total observing times on target were 10185 s for the 
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complete 2.38 - 45 /mi grating scan and 3962 s for the three Fabry- Perot scans. 

3. Results 

The reduction of the grating spectrum was carried out using the SWS Interactive Analysis 
system, applying calibration files from version 7.0 of the ISO pipeline software. The Fabry-Perot 
spectra were initially reduced using version 6.22 of the ISO pipeline software, and the ISAP 
software package^] was then used to remove bad data points and to co-add the individual spectral 
scans. 

The 29.5 — 45 fim grating scan of VY CMa, shown in Figure 1, reveals a strikingly rich 
spectrum of water rotational transitions. At least 41 features are readily identifiable as water, 
corresponding to a density of more than 2 such features per micron of wavelength coverage. Line 
identifications and flux estimates are given in Table 1. Here we have been quite conservative in 
identifying spectral features with water transitions, including only those features for which the 
identification seems most secure. 

The flux calibration is believed accurate to ~ 30%, and the typical la statistical error in the 
flux is ~ several x 10 -20 Wcm -2 . Multiple scans of some portions of the spectrum strongly suggest 
that all features apparent to the eye are real. However, many features, whether in absorption or 
emission, are still not reliably identified (although many of them may, in fact, be water vapor 
lines). As a result, the proper identification of a baseline, on top of which individual features 
are observed, is quite uncertain. For many water lines, this difficulty substantially increases the 
uncertainties in the line flux estimates given in Table 1. 

In addition to transitions within the ground vibrational state, the line list in Table 1 includes 
rotational transitions within the (010) excited vibrational state. Two radio wavelength transitions 
within the (010) state had previously been detected in VY CMa by Menten & Melnick (1989). 

The only molecule other than water for which a feature has been identified in this spectral 
region is OH, the primary photodissociation product of water, which we detected by means of the 
2 n 1/2 (J = 5/2) n 3/2 (J = 3/2) absorption feature at 34.6 /mi. 

Several broad spectral features have been identified in the 20 - 45 /mi spectra of other evolved 
stars (e.g. Waters et al. 1996), some of which have been tentatively attributed to crystalline 
silicates. One broad spectral feature is prominent in the 29.5 - 45 /mi VY CMa spectrum (Figure 
1). It is centered at wavelength 30.4 ± 0.1/mi, with a width ~ 0.5 /im, and a peak flux lying ~ 3% 
above the continuum. A feature with very similar properties has been detected in each of the 
sources NML Cyg, AFGL 4106, HD 179821 and NGC 6302 (Waters et al. 1996): its identification 
is unknown at present. 



lr The ISO Spectral Analysis Package (ISAP) is a joint development by the LWS and SWS Instrument Teams and 
Data Centers. Contributing institutes are CESR, IAS, IPAC, MPE, RAL and SRON. 
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Figure 2 shows the Fabry-Perot spectra obtained toward VY CMa. The corresponding line 
fluxes are listed in Table 2: once again, the likely errors in the flux calibration are believed to be 
~ 30%. 

4. Discussion 

4.1. Grating spectrum 

The grating spectrum of VY CMa shown in Figure 1 demonstrates dramatically the ubiquitous 
nature of water vapor emission lines in the far-infrared spectral region. Adopting 1500 pc as the 
most credible estimate of the distance to VY CMa (Lada Sz Reid 1978), we find that the identified 
water features listed in Table 1 carry a total luminosity of 25 Lq. Detailed modeling of the 
water line fluxes, which we defer to a later publication, will be complicated significantly by large 
departures from spherical symmetry. In particular, recent Hubble Space Telescope observations 
of VY CMa (Kastner & Weintraub 1998) vividly demonstrate the highly anisotropic nature 
of the surrounding outflow. These observations suggest that a highly flattened circumstellar 
envelope obscures the photosphere of VY CMa from direct view along our line of sight, but that 
photospheric radiation does leak out along other directions, giving rise to an elongated reflection 
nebula. Clearly, this complex geometry is very different from the idealized spherical geometry 
that has been assumed previously in modeling water emissions from circumstellar outflows (e.g. 
Deguchi & Rieu 1990, Chen & Neufeld 1995). We expect that future modeling efforts will be 
usefully constrained by additional observations of the water vibrational bands near 2.7 and 6.2 fim 
and of rotational absorption lines in the 16 - 19.5 fim spectral region, the reduction and analysis 
of which are underway. 

The strength of the 34.6 /im OH absorption feature is of relevance to models for the radiative 
pumping of 1612 MHz OH masers by 34.6 fim infrared radiation (e.g. Elitzur, Goldreich & 
Scoville 1976). The OH absorption that we observed presumably arises in a shell surrounding the 
source where the outflowing water vapor is photodissociated by the interstellar radiation field. 
Assuming spherical symmetry, we find that the total rate of absorption in the OH 34.6 fim feature 
is 1.0 x 10 47 photons per second (for an assumed distance to the source of 1500 pc). The average 
maser photon emission rate is 4.3 x 10 45 s" 1 (Harvey et al. 1974), with a typical time variability 
;$ 20%. Disregarding the complications introduced by anisotropy, we therefore find that ~ 23 
photons are absorbed in the 34.6 ^m 2 U 1/2 (J = 5/2) U 3/2 (J = 3/2) line for every 1612 MHz 
maser photon emitted. The ratio of 34.6 ^m photons absorbed to 1612 MHz photons emitted is 
somewhat higher than the value of 14 reported by Sylvester et al. (1997) for the evolved star of 
IRC+ 10420. 
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4.2. Fabry-Perot spectra of VY CMa 

The Fabry-Perot spectra of VY CMa shown in Figure 2 are of sufficiently high spectral 
resolution to yield useful information about the line profile. As far as we are aware, these are the 
first spectra in which the line profiles of thermal water emissions from a circumstellar outflow have 
been resolved spectrally. All three observed lines exhibit the classic "P Cygni" profiles associated 
with line formation in an outflow. We find that the observed line profiles are consistent with 
the line emission expected from a spherical outflow at constant velocity, t> ou t> accompanied by 
continuum emission that suffers blueshifted absorption in the approaching side of the outflow. 
Figure 2 shows that an outflow velocity of t> out = 25kms~ 1 yields an acceptable fit for all three 
lines: here the dashed line represents the symmetric, parabolic emission line profile expected for 
optically-thick line emission from a constant velocity outflow; the dotted line represents continuum 
radiation with an optically-thick absorption feature blueshifted by velocity v ou t relative to the the 
emission line; and the solid line represents the sum of those components after convolution with the 
instrumental profile. 

The integrated fluxes in the emission and absorption features, the equivalent widths of the 
absorption features, and the emission line velocity relative to the LSR, ve, are all tabulated in 
Table 2. The values of ve are in acceptable agreement with each other and with the stellar LSR 
velocity of ~ 18kms _1 inferred for this source (Reid & Dickenson 1976). The value of 25kms~ 1 
obtained for f out is close to - but somewhat smaller than - the estimates of ~ 37kms _1 and 
~ 32kms _1 inferred previously for the terminal outflow velocity from observations of the thermal 
J = 2 - 1 line of SiO (Reid & Dickenson 1976) and the 1612 MHz OH maser line (Reid & 
Muhleman 1978). The observed spectra therefore provide an important observational confirmation 
of theoretical predictions (e.g. Deguchi & Rieu 1990, Chen & Neufeld 1995) that the far-infrared 
water line emission detected from circumstellar outflows is generated largely within a region where 
the outflow has already reached - or nearly reached - its terminal velocity. The absorption line 
equivalent width is considerably larger for the 432 — 3o3 line at 40.6909 fim than for the other two 
lines, a difference that presumably results from a large optical depth in this transition which has a 
lower state of energy, Ei/k, only 197 K. 
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and NAG5-3542 and NASA contract NAS5-30702 to G.J.M. M.H. also thanks the Alexander 
von Humboldt foundation and the Max Planck Institute for Radioastronomy for their hospitality 
during a three-month visit to the MPIfR in Bonn. 
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TABLE 1 

Lines identified in the grating scan of VY CMa 



Line Wavelength Upper state Flux 

Energy (E/k) 

(pm) (K) (1Q-18 Wcm- 2 ) 



H 2 7 6i - 652 


30.526 


1750 


0.91 


H2O 762 — 651 


30.529 


1750 




H2O 44i — 3i2 


31.772 


702 


0.65 


H 2 643 - 616 


32.313 


1089 


0.88 


H 2 7 5 2 - 643 


32.991 


1525 


2.22 


H 2 6 6 o - 5 5 i 


33.005 


1504 




H 2 661 - 5so 


33.005 


1504 




H 2 IO47 - 936 


33.510 


2275 


0.90 


H2O 5so — 523 


33.833 


1068 


0.82 


H2O 946 — 835 


34.396 


1929 


0.68 


H2O 734 — 625 


34.549 


1212 


0.58 


OH J = 5/2- «— 3/2+ b 


34.603 


416 


-1.49 


OH J = 5/2+ <— 3/2- b 


34.629 


415 


-0.61 


H2O 743 — 634 


35.429 


1340 


1.85 


H 2 5 33 - 4 04 


35.471 


725 


1.52 


H2O 845 - 734 


35.669 


1615 


0.94 


H2O 651 - 542 


35.904 


1279 


0.60 


H 2 652 - 5 4 i 


35.938 


1279 


1.21 


H 2 7 5 2 - 725 


36.046 


1525 


0.97 


H 2 V2 5so - 44i 


36.161 


3462 


0.74 


H 2 u 2 5 5 i - 4 40 


36.163 


3462 




H2O 624 — 5l5 


36.212 


867 


0.56 


H2O 123 i io — II29 


36.786 


2824 


0.40 


H2O 744 - 633 


37.566 


1335 


0.90 


H2O 441 - 4l4 


37.984 


702 


0.80 


H 2 II39 - 10 2 8 


38.895 


2439 


0.20 


H2O 7 3 4 - 7()7 


39.045 


1212 


0.67 


H 2 5 50 - 4 41 


39.375 


1068 


2.31 


H2O 5si - 440 


39.380 


1068 




H 2 642 - 5 3 3 


39.399 


1090 




H 2 12 2 ,n - Hi, 10 


40.016 


2554 


0.40 


H 2 844 - 817 


40.179 


1628 


0.77 


H 2 13l,l 3 - 12o,12 


40.188 


2600 




H 2 13 ,13 - 12l,l2 


40.189 


2600 




H2O 643 — 532 


40.337 


1089 


1.32 


H 2 u 2 541 - 4 32 


40.478 


3240 


0.90 


H2U U2 042 — 431 


40.687 


3240 


2.21 


H2O 4 32 - 3 03 


40.691 


550 


a 


H2O 633 - 524 


40.760 


952 


1.17 


H 2 1129 - 1038 


40.894 


2433 


0.54 


H 2 10 38 - 9 27 


40.948 


2081 


0.63 


H2O 132,12 — 13l,13 


42.427 


2939 


0.29 


H 2 13l,12 - 13 ,13 


42.438 


2938 


a 


H 2 9 37 " 8 26 


42.859 


1750 


0.79 


H2O V2 523 - 4l4 


43.035 


2955 


0.65 


H2O 1138 - IO47 


43.125 


2609 


0.50 


H 2 111, 10 - IO29 


43.250 


2194 


0.78 


H 2 12i,i2 - 110,11 


43.339 


2241 


0.24 


H2O 12o,l2 — 111, 11 


43.341 


2241 


a 


H 2 U 2 73 5 - 624 ? 


43.713 


3510 


0.49 


H 2 541 - 4 32 


43.893 


878 


1.35 


H 2 7 43 - 7i 6 


44.048 


1340 


1.15 


H 2 5 42 - 4 3 i 


44.195 


878 


1.14 


H2O 836 - 725 


44.702 


1448 


1.19 



a Line blended with previous feature. Quoted flux represents a sum over all blended lines. 
b Cross-ladder 2 n 1 / 2 4—2 n 3 / 2 transitions. 



-9- 



TABLE 2 

Results of Fabry-Perot observations of VY CMa 



Line 




Rest 


Upper state 


Line Fluxes 


Absorption line 


V E ° 






Wavelength b 


Energy (E/k) 


Emission Absorption 


Equiv. Width 












F E F a 


W A 








W 


(K) 


(10- 18 Wcm- 2 ) 


(km s- 1 ) 


(km s- 1 ) 


H 2 7 25 


- 6l6 


29.8367 


1126 


2.35 -1.94 


11 


22 


H 2 4 4 i 


-3l2 


31.7721 


702 


2.41 -1.71 


11 


20 


H 2 4 32 


— 3o3 


40.6909 


550 


3.69 -2.20 


32 


18 



a Central velocity of the emission feature, relative to the LSR 
b From Toth (1991) 
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Fig. 1. — 29.5 - 45 fim SWS grating scan of VY CMa. With the exception of the OH absorption 
feature at 34.6 /xm, all the identified features are water lines. 
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Fig. 2. — SWS Fabry-Perot spectra of three pure rotational lines of water observed towards VY 
CMa: the 7 25 - 6i 6 line at 29.8367/mi, the 4 4i - 3 i2 line at 31.7721 /xm, and the 4 32 - 3 3 
line at 40.6909 fim. Solid curves show the fits described in the text, obtained by convolving the 
instrumental profile with the sum of the two components represented by the dotted and dashed 
lines. 



